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I.  STATEMENT  OF  THE  PROBLEM  AND  APPROACH 

Models  of  the  combustion  of  nitrate  ester  and  nitramine  propellants,  especially 
those  for  dark  zone  structures,  depend  critically  on  chemical  kinetic  input  data  for 
individual  reactions  at  the  temperatures  of  interest.  With  the  ARO  support  we  have 
obtained  such  data  for  sensitive  reactions,  in  isolation  of  other  reactions,  at  realistic 
temperatures  for  dark  zones.  One  of  our  HTP  (high-temperature  photochemistry) 
reactors,  which  uniquely  allow  such  studies  to  be  made,  was  used.  Furthermore,  an 
extensive  book  chapter  was  prepared.1  It  (i)  reviews  the  current  knowledge  of  the 
gaseous  phase  chemistry  of  propellant  combustion,  (ii)  discusses  the  rate  coefficients  and 
mechanisms  of  the  important  individual  reactions,  and  (iii)  describes  the  techniques  by 
which  this  kinetic  input  information  is  obtained  for  the  relevant  temperatures. 

II.  SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 
A.  Influence  of  chemical  kinetics  on  propellant  combustion  models. 

The  book  chapter1  extensively  discusses  the  influence  of  chemical  kinetic  data  on 
propellant  burning  models.  Included  in  this  are  several  reactions  which,  if  included  with 
the  prior  literature  rate  coefficients,  suggested  behavior  fastly  different  from  what  is 
observed.  This  is  illustrated  in  Fig.  1  for  the  CO?  production  from  Hj/CO/NO/Ar 
mixtures  at  1000  K.  H2,  CO,  and  NO  are  all  major  components  of  nitrate  ester  propellant 
dark  zones,  DZ.  Fig.  la  shows  the  experimental  points  obtained  by  Diau  et  al.  for  such 
mixtures  at  1000  K.  Also  shown  in  the  figure  is  the  calculated  curve  for  the  CO2 
production  using  the  ARL  (Dr.  W.R.  Anderson)  baseline  DZ  model.  Very  good 
agreement  is  obtained.  However,  if  prior  literature  data  for  each  of  three  reactions  and, 
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because  of  microscopic  reversibility,  their  reverse  processes,  are  included  in  the  models, 
the  results  would  be  at  variance  with  the  observations. 

Fig.  lb  shows  the  result  when  the  first  of  these  reactions 

N  +  C02  =  N0  +  C0  (1) 

is  added  to  the  baseline  mechanism,  by  using  the  experimental  forward  rate  coefficients 
k(T)  from  Ref.  3  and  thermodynamics  to  compute  the  rate  coefficient  of  the  assumed 
reverse  reaction.  The  calculated  C02  production  rate  now  would  be  two  orders  of 
magnitude  too  fast.  If  instead 

NH  +  H20  =  HNO  +  H2  (2) 

is  added,  using  the  forward  k(T)  of  Ref.  4  and  the  thermochemistry,  the  predicted  rate  is 
about  three  orders  too  fast,  Fig.  lc.  And  if  the  forward  k(T)  for4 

NH  +  C02  =  HNO  +  CO  (3) 

were  similarly  added,  the  rate  becomes  four  orders  of  magnitude  too  fast.  The  result  of 
using  these  literature  values3'4  in  DZ  models  is  that  the  predicted  chemical  delay  time 
would  be  much  too  short  and  the  DZ  length  would  be  reduced  such  that  the  DZ  would 
disappear.  It  therefore  appeared  that  either,  case  a,  the  reported  rate  coefficients  for 
reactions  (2-4)  are  too  large,  or  case  b,  the  products  are  not  those  suggested  by  these 
reaction  equations.  To  check  these  possibilities  the  chemical  kinetic  studies  summarized 
in  Section  II  B  were  performed. 

B.  Chemical  Kinetic  Studies 

Both  experimental  and  theoretical  investigations  have  been  made.  The  study  of 
reaction  (1)  was  done  under  a  preceeding  ARO  grant.  No  reaction  was  observable, 
kj(285-  1 142K)  <  5  x  10’16  cnrmolecule'V1,5  allowing  the  reaction  to  be  excluded  from 
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the  models.  The  results  of  the  Avramenko  study3  were  probably  due  to  the  authors’ 
erroneously  attributing  the  faster  N-atom  recombination,  resulting  from  the  presence  of 
CO2  as  an  inert  third-body,  to  a  reaction  with  this  species.  As  a  check  on  the  HTP 
method  used  in  our  work,  another  N-atom  reaction  was  studied.  For  this  reaction,  N  +  O2 
=  NO  +  O,  finite  reaction  rate  coefficients  were  measured.  In  the  lower  temperature 
range  of  overlap  with  earlier  works  these  showed  excellent  agreement.  A  subsequent  ab 
initio  study  of  reaction  (1)  at  ARL  showed  that  the  barrier  in  the  lowest  possible  path  is 
so  large  that  no  reaction  could  have  been  observed.6  Thus,  reaction  (1)  is  an  example  of 
case  a  above. 

Reactions  (2)  and  (3)  have  been  found  in  the  present  work7  to  be  examples  of  case 
b.  In  the  Rohrig  and  Wagner  study4  NH  was  produced  in  shock  waves  by  HN3 
dissociation.  To  verify  their  rate  coefficients,  we  first  studied  the  other  reaction  reported 
in  their  work 

NH  +  H2  =  NH2+H  (4) 

In  the  HTP  reactor  the  NH  was  produced  by  the  2-photon  193  nm  photolysis  of  NH3.  We 
obtained  k(833- 1432  K)=3.5  x  10'uexp(-7758  K/T)  cm3 molecule'V1.  This  is  in  such 
excellent  agreement  with  their  work  from  1280  to  1685  K  that  the  results  could  be 
combined  to  yield  k(833- 1 685  K)=4.4  x  10~nexp(-8142  K/T)  cm3 molecule'V1.  This 
made  it  highly  unlikely  that  there  could  be  a  serious  error  in  the  k2  and  k3  measurements 
of  Ref.  4,  thus  indicating  that  the  products  are  probably  not  HNO  +  H2  and  HNO  +  CO, 
respectively.  Direct  measurements  of  k2  and  k3  by  the  HTP  method  did  not  work  out  as 
the  photolysis  produced  O  atoms,  which  react  rapidly  with  the  NH.  It  was  therefore 
decided  to  study  one  of  these  reactions,  i.e.  NH  +  C02,  by  ab  initio  means. 
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Details  of  the  theoretical  calculations  are  also  discussed  in  Ref.  7.  The  derived 


potential  energy  diagram  is  shown  in  Fig.  2.  The  dotted  line  gives  the  triplet  surface  path 
for  the  ground  state  triplet  3NH.  Production  of  HNO  +  CO  by  this  route  involves  an 
energy  barrier  of  226  kJ  mol1 .  Alternately,  a  spin-flip  could  occur  to  a  singlet  surface 
path.  This  channel,  involving  the  electronically  excited  *NH  molecule,  is  shown  by  the 
solid  line.  This  still  requires  passage  over  a  126  kJ  mol  1  barrier  to  produce  HNO  +  CO. 
As  the  measured4  activation  energy  is  only  60  kJ  mol'1  it  may  be  concluded  that  HNO  + 
CO  are  not  the  reaction  products. 

In  an  attempt  to  establish  what  the  actual  products  of  the  NH+CO2  reaction  are, 
ab  initio  enthalpies  for  the  formation  of  several  adducts  were  calculated.  The  most  stable 
of  these  HONCO  is  exothermic  by  only  33  kJ  mol'1,  which  is  insufficiently  stable  to  be 
the  final  reaction  product.  It,  or  another  adduct,  likely  will  react  further  with  other 
species  present.  Reaction  (2)  likewise  probably  involves  initial  adduct  fonnation.  This 
was  not  further  investigated  here.  It  should  be  mentioned  that  H3NO  and  H2NOH  have 
been  suggested4  as  such  energetically  accessible  products  of  that  reaction. 

The  chemistry  problem  of  the  identity  of  the  products  of  reactions  (2)  and  (3)  thus 
remains.  However,  it  can  be  considered  as  settled  that  reactions  (1-3),  with  their  reverse 
processes,  should  not  be  included  in  DZ  models. 
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Fig.  1.  Comparison  of  experimental  and  modeling  results  for  an  F^/CO/NO/Ar  mixture 

at  1000  K.  Note  the  different  timescales  in  (a)-(d). 

(a)  Experimental  results  (points)  and  modeling  result  using  the  baseline  mechanism 
(curve). 

(b)  Modeling  result  when  N  +  CO2  =  NO  +  CO  (1)  is  added  to  the  baseline 
mechanism,  using  the  experimental  forward  k(T)  from  Ref.  3  and 
thennodynamics  to  compute  the  rate  of  the  assumed  reverse  reaction. 

(c)  Modeling  result  when  NH  +  H2O  =  HNO  +  Ff>  (2)  is  added  to  the  baseline 
mechanism,  using  experimental  forward  k(T)  from  Ref.  4  and  thermodynamics  to 
compute  the  rate  of  the  assumed  reverse  reaction. 

(d)  Modeling  result  when  NH  +  CO2  =  HNO  +  CO  (3)  is  added  to  the  baseline 
mechanism,  using  the  experimental  forward  k(T)  from  Ref.  4  and 
thennodynamics  to  compute  the  rate  of  the  assumed  reverse  reaction. 
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Enthalpy  relative  to  NH  +  C02  at  298  K,  kJ  mol'1 


Fig.  2.  Potential  energy  diagram  for  NH+CO2  showing  paths  to  HNO+CO.  The  solid 
line  shows  the  singlet  path,  the  dashed  line  the  triplet  path. 
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